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ABSTRACT 

We present accurate absolute astrometry of 6.7 GHz methanol masers de- 
tected in the Arecibo Methanol Maser Galactic Plane Survey using MERLIN 
and the Expanded Very Large Array (EVLA). We estimate the absolute astrom- 
etry to be accurate to better than 15 and 80 milliarcseconds for the MERLIN 
and EVLA observations respectively. We also derive the morphologies of the 
maser emission distributions for sources stronger than ~ 1 Jy. The median 
spatial extent along the major axis of the regions showing maser emission is 
~ 775 AU. We find a majority of methanol maser morphologies to be complex 
with some sources previously determined to have regular morphologies in fact 
being embedded within larger structures. This suggests that some maser spots 
do not have a compact core, which leads them being resolved in high angular 
resolution observations. This also casts doubt on interpretations of the origin of 
methanol maser emission solely based on source morphologies. We also investi- 
gate the association of methanol masers with mid-infrared emission and find very 
close correspondence between methanol masers and 24 /im point sources. This 
adds further credence to theoretical models that predict methanol masers to be 
pumped by warm dust emission and firmly reinforces the finding that Class II 
methanol masers are unambiguous tracers of embedded high-mass protostars. 

Subject headings: masers — instrumentation: high angular resolution — astrom- 
etry — stars: formation 
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Introduction 



The 6668.519 MHz Jk = 5n — 61 transition is the strongest and most widespread 



Class II methanol maser line (jMentenlll99ll ). There is strong evidence that 6.7 GHz methanol 
masers trace very early evolutionary phases of high-mass star formation. For example, 
multi- wavelength studies of 6.7 GHz me thanol masers have found their environment to be 
characteristic of massive protoclusters (IMinier et al.l 120051 ). Moreover, recent work mod- 
eling the spectral energy distributions of the masers' host sources have found them to be 
consistent with that of massive young stellar objects (MYSOs) undergoing rapid accretion 
(IPandian et aPboiol ). 



However, the question where the maser emission ori ginates relative to the MYSO is 



still t he subject of considerable debate. Early studies (e.g. iNorris et al.lll993l : iPhillips et al. 
19981 ) found a number of methanol masers to have linear or arched morphologies with mono- 
tonic velocity gradients across the spatial features. These were interpreted to be indicative 
of the masers originating in an edge-on circumstellar disk. However, the minimum mass 
for the central o bject derived from the disk hypothesis was found to be unrealistically high 
in some sources (jWalsh et al.l 119981 ). In contrast, assuming that the linear maser distribu- 
tions trace the full extent of such a disk, the enclosed mass was found to be much less 
than 1 Mq in a number of sources ( iMinier et al.l I2OOOI ) . This led to suggestions that the 
masers only trace partial disks a round massive stars, or are associated with outflows or 
shock fronts. iDodson et al.l (|2004[ ) suggested that methanol masers originated behind a pla- 
nar shock propagating through a star forming core in order to explain the observation of 
velocity gradients within individual spot clusters that were perpendicular to the main large 
scale velocity gradient across the linear morphology. 

High spatial resolution mid-infrared studies suggest that many methanol masers, includ- 
ing some with linear morphologies, are associated with outflows rather than being located 
in disks. While G 35.20-0.74 is th e best example of masers occurring along the walls of 
an outflow cavity (jPe Buizerl |2006| ) , seve ral other masers were found to have their linear 



such as NGC 7538 IRSl flPestalozzi et al. 



morphologies aligned par allel to outflows ( De Buizerl 120031 ). However, there are also sources 



20091 



\ where there is strong evidence that masers 



originate in a disk. Furthermore, in the recent study of the proper motions of 12.2 GHz 
methanol masers (whic h are the second strongest Class II methanol masers) in W3(0H), 
Moscadelli et al.l ( 120101 ) find that the linear distribution of the masers can be well fit with 
a fiat rotating disk that is seen almost edge-on. Since there is very good spatial correlation 
between the locations of 6.7 and 12.2 GHz methanol masers, this can be taken as evidence 
for a disk origin of some 6.7 GHz methanol masers. 



Recent high resolution imaging of 6.7 GHz methanol masers with the European VLSI 
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network (EVN) revealed a number of masers to have ring shaped morphologies ( iBartkiewicz et al, 
20091 ). The ring morphologies were modeled as an inclined disk or torus around a MYSO 
having expanding or infalling kinematics. In addition, a number of masers were found to be 
associated with sources that had excess emission mea sured in the 4.5 /im IRAC band with the 
Spitzer space telescope (see ICyganowski et al.ll2009l for a detailed discussion of these sources 
called Extended Green Objects or EGOs). Taken together, this was interpreted as evidence 
for the origin of methanol maser emission in the interface region between a disk/torus and 
an outflow. 

One of the potential problems in very high angular res olution studies is th at a significant 
fraction of the maser flux is resolved out. For example, iMinier et al.l (|2002| ) found that a 
number of maser spots had a compact core and a more diffuse halo, while some spectral 
features had no detectable compact core. Hence, lower angular resolution observations are 
required to determine the full extent and morphology of maser emission. It is also of interest 
to carry out this work towards a large and homogeneous sample from a blind survey. 

The most sensitive blind survey to date for 6 .7 GHz methanol rn asers is the Arecibo 
Methanol maser Galactic Plane Survey (AMGPS; iPandian et al.l 120071 . hereafter Paper I). 
AMGPS detected 86 6.7 GHz methanol masers, of which 48 were new detections. Taking 
into account the pointing accuracy of the Arecibo telescope, th e uncertainty in the position 
of the AMGPS masers is ~ 18" at the 95% confidence level (IPandian fc Goldsmith! l2007l . 
Paper II hereafter). Paper II also presented an analysis of the properties of the masers 
including association with mid-infrared sources, a number of methanol masers were found to 
not have mid-infrar ed counterparts in t he Galactic plane sur vey using the Midcourse Source 
Experiment (MSX: iPrice et al.ll200ll ). IPandian et al.l (120091 Paper HI hereafter) presented 
the systemic velocities and distances to the sources (all the distances used in this paper are 
taken from Paper HI) leading to the determination of the luminosity function of 6.7 GHz 
methanol masers. In this paper, we present accurate astrometry for the AMGPS sources 
usmg MERLI^Q and the Expanded Very Large Arraj]^ (EVLA). We also derive morphologies 
of the maser emission for the strong sources (peak flux densities > 1 Jy), and examine the 
association of methanol masers with mid-infrared sources found in recent sensitive surveys 
using the Spitzer space telescope. 



^ Based on observations made with MERLIN, a National Facility operated by the University of Manchester 
at Jodrell Bank Observatory on behalf of STFC. 

^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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2. Observations and data reduction 

2.1. MERLIN observations 

The MERLIN observations of the AMGPS masers (W51 main and G41.87-0.10 ex- 
cluded) were carried out between 2007 February and June using six antennas except for 
observations between March 5 and March 11 which employed only five antennas. Two corre- 
lator configurations were used for the observations. In the wide-band mode, two polarizations 
were observed with a total bandwidth of 16 MHz split into 32 channels with a channel width 
of 0.5 MHz. In the narrow-band mode, two polarizations were observed with a bandwidth 
of 1 MHz split into 256 channels yielding a channel spacing of 3.9 kHz. The primary and 
bandpass calibrator, 3C84 was observed in both wide and narrow-band modes. The phase 
calibrators (see Table [T]) were only observed in wide-band mode to achieve adequate signal 
to noise ratio, while the target sources were only observed in narrow-band mode. Phase 
referencing was used with each scan on a target source being preceded and followed by a 
scan on a phase calibrator. To achieve good uv-coverage, the target observing sequence was 
cycled through a group of sources with similar systemic velocities. No Doppler tracking was 
used, and each group of targets was observed at a fixed frequency. This observing sequence 
was repeated several times over the course of each observing run which typically lasted for 
~ 10 hours. The final integration time was roughly 45 minutes per source, with the resulting 
la channel noise being ~ 45 mJy beam~^. At the adopted rest frequency of 6668.519 MHz 
for the maser transition, the velocity resolution in the narrow-band mode is 0.18 km s~^. 



The data were initially reduced using local MERLIN software (jPiamond et al.ll2003l ). 
After initial editing and flux calibration on 3C84, the data were converted into FITS for- 
mat with subsequent processing being carried out using the Astronomical Image Processing 
System (AIPS) package of NRAO. The flux calibration was carried out assuming the flux 
density of 3C84 to be 14.5 Jy in 2007 February, 15.0 Jy in 2007 March, 15.5 Jy in 2007 
April, 16.0 Jy in 2007 May, and 16.5 Jy in 2007 June. The phases of 3C84 were calibrated 
in both the wide-band and narrow-band data, with the derived phase offset between the two 
data sets being used to transfer the calibration from the phase reference source (which is 
only observed in the wide-band mode) to the target sources (which are only observed in the 
narrow-band mode). After calibration, the task "CVEL" was used to shift the spectra to 
the systemic velocity with respect to the local standard of rest (LSR) and correct for the 
effects of Earth's motion. The target spectra were then inspected followed by imaging the 
brightest channel to obtain the target astrometry. The typical beamsize was 60 x 35 mas 
with a position angle of 20°. 



As indicated in Sect. 2.3.1 of 



Bartkiewicz et al, 



(120091 ) ■ the astrometric accuracy is 
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limited by four factors - the accuracy of the phase cahbrator position, accuracy in the 
positions of the individual antennas, accuracy in the transfer of phase solutions from the 
phase calibrator and targets (on account of the atmosphere), and the precision with which 
the emission location can be determined given the beamsize. The phase calibrator positions 
indicated in Table [1] have position uncertainties between 0.4 and 4.4 milliarcseconds (mas). 
The antenna positions have uncertainties of 1-2 cm, which results in a position uncertainty of 
~ 10 mas. To determine the uncertainties from the phase transfer from the phase calibrator 
to the targets, we inspected the phase change over a timescale corresponding to the separation 
between the calibrator and target (e.g. an angular separation of 3° corresponds to a 12 min 
timescale). We estimate positional uncertainties to be less than 10 mas, with median values of 
6 mas. The last uncertainty from formal fitting errors is usually <C 1 mas, and can usually be 
ignored. Combining the various uncertainties, we estimate the absolute astrometric accuracy 
to be better than 15 mas. 



2.2. EVLA observations 

Among the 81 AMGPS masers targeted in the MERLIN observations, 30 sources were 
either not detected, or had bad data. After verification of the peak flux densities of these 
sources with the 100 m Effelsberg telescope to check for variability, we observed 27 sources 
in 24 fields using the EVLA. The observations were carried out on 2008 October 16 and 18 in 
the A configuration. A single polarization (RR) was recorded using the old VLA correlator 
with a bandwidth of 1.56 MHz and 512 spectral channels. Since the lower 0.5 MHz of the 
bandpass was unusable due to aliasing, this setup was required to achieve the highest possible 
spectral resolution along with the required velocity coverage. The resulting channel spacing 
was 3 kHz giving a velocity resolution of 0.14 km s~^ at the rest frequency of the maser 
transition. 

Each target source was observed in a single 10 minute or two 15 minute scans depending 
on its single dish flux density, with the resulting Icr channel noise ranging from 20 to 35 mJy 
beam~^. The sky frequencies were computed for each individual scan and placed within 
the upper two-thirds of the 1.56 MHz bandwidth to ensure that no maser component would 
fall in the lower 0.5 MHz of the bandpass. The source J1331+305 (3C286) was used to 
set the absolute flux density scale, and the sources J1751+096 and J1925+211 as bandpass 
calibrators. The phase calibrators, used to calibrate the complex gains were J1851+005, 
J1856+061 and J1922+155 (Table [2]), with each calibrator being used for a subset of the 
target sources closest in angular distance. 

The editing, calibration, and processing of the data were carried out using AlPS. After 
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the transfer of the complex gain cahbration solutions to the target sources, the strongest 
emission channel in each source was imaged and its position measured followed by self- 
calibration (when sufficient signal to noise was available) in both phase and amplitude in 
an iterative cycle. In several cases, multi-field imaging was necessary to account for strong 
6.7 GHz methanol maser emission from other sources that fell within the primary beam of 
the 25 m EVLA antennas. In some weak sources, multiple channels were smoothed in order 
to detect the source. 

The phase calibrators used in the EVLA observations have a position accuracy better 
than 3 mas (Table [2]). The EVLA antenna positions are known to an accuracy of ~ 1 cm 
implying a position uncertainty of 55 mas. The uncertainty in the astrometry resulting 
from the transfer of phase solutions from the phase calibrator to the target (obtained from 
examining the raw phases of the phase calibrator) is estimated to be less than 30 mas. In 
addition, some sources are only detected with moderate to poor signal to noise ratio, and 
consequently the formal fitting errors in determining the position of peak emission are as 
high as ~ 45 mas for these sources. Thus, we estimate the absolute astrometry to be accurate 
to between 65 and 80 mas. 



2.3. 24 fiui counterparts 

To examine the relation between 6.7 GHz methanol masers and mid- infrared emission 



we looked for 24 fim counterparts in the Spitzer MIPSGAL survey (jCarey et al.ll2009l ). Since 
the 24 (Um point source catalog has not been released, we used the mosaiced images corrected 
for artifacts for this purpose. The strong interstellar dust emission in the Galactic plane leads 
to complex backgrounds in the images, which can significantly affect the measurement of the 
positions of point sources. Hence, we used a 11 x 11 pixel median filter to subtract the 
background, as suggested in the MIPSGAL data reduction cookbooljf]. The use of such 
aggressive filtering reduces the source fiux to some extent. While this can be accounted for, 
we are primarily interested in the positions of the point sources in this work. Hence, we do 
not derive fiux densities of the 24 fim counterparts. 

An additional complication arises from many regions being extremely bright and satu- 
rated in the MIPSGAL data. While some regions are completely blanked and other sources 
have the entire core of the point spread function (PSF) blanked, there are some sources 
for which only the central few pixels are blanked. Although the latter are rejected in the 
point source detection algorithm of MOPEX, there is sufficient information to measure the 



■^http: / /ssc. spitzer. caltech.edu/dataanalysistools /cookbook / 30/#_Toc272753806 
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positions of these sources. Hence, we measured the positions of all 24 fj,m counterparts by 
Gaussian fitting of the PSF cores using the routine "JMFIT" in AIPS. Although the full 
width at half maximum (FWHM) of the MIPS PSF at 24 ^m is 5.9", the high signal to 
noise ratio of the counterparts implies that the formal uncertainties in the Gaussian fits are 
small, and that systematic errors dominate the uncertainties in the positions derived for the 
MIPSGAL counterparts. 



3. Results 



We were able to measure the absolute positions of 57 sources u sing MERLIN (Table [31 
positions of masers in W51 main were obtained using the data of IXu et al.l |2009| ) and an 
additional 25 sources using the EVLA (Table H]). We also identified two masers in the 
MERLIN d ata that were not in t he original AMGPS catalog (one of these sources was 
detected by lCyganowski et al.ll2009l ). Both sources are located close to other sources, leading 
them to be either missed or unresolved in the original survey. Thus, this work has determined 
the absolute positions to 82 out of 88 AMGPS sources. 

Point source counterparts at 24 /zm were found for 60 out of 82 sources. 18 sources 
were either fully saturated or in the middle of a blanked region. Another two sources was 
found to be in the middle of extended emission, while no counterparts were found for two 
sources. The positions of the 24 /zm point source counterparts are tabulated in Table [5l and 
comments about selected sources are included in Sect. 13.21 



3.1. Relative astrometry and spot morphology 

Although the absolute astrometric accuracy is limited by the factors described in Sect. 
12.11 and 12.21 the relative astrometry between the emission distributions in different velocity 
channels can be determined to much higher precision. To this end, we imaged each source 
over the full velocity range of its emission, and measured the positions of the maser spots 
relative to the channel with peak emission. To improve the dynamic range in the data cubes, 
we self-calibrated data of targets that were sufficiently strong (peak flux density > 1 Jy). 
The measurement of spot positions in the data cube was automated using the AIPS task 
"SAD" . To obtain reliable fits, we imposed a 6a threshold on both the peak and integrated 
fiux density to qualify as a maser spot. 

Figs. [T]and[2]show the results of this work for the MERLIN and EVLA data respectively. 
For each source, the top panel shows the integrated spectrum, while the bottom panel shows 
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the morphology of the maser spot emission as a function of both angular and physical 
distance. The coordinates listed in Tables [3] and H] serve as reference coordinates for the spot 
maps. The LSR velocities of the maser spots are indicated next to the spots themselves. 
When a group of spots form a morphology with a monotonic velocity gradient, only the initial 
and final velocities are indicated at the ends of the feature. In cases where there is a compact 
cluster of spots, or where the spots form a morphology without a monotonic velocity gradient, 
the velocity range of the cluster is indicated next to the feature. Some sources observed with 
MERLIN are detected in three spots or fewer, from which no morphological information can 
be deduced. For these cases, only the integrated spectra are shown in Fig. [31 The EVLA data 
have much poorer resolution compared to the MERLIN data, and thus spot morphologies can 
be deduced only for the strong sources, or for sources that show a significant spatial extent 
of spot emission. Hence, for all EVLA sources where the maser spots are concentrated in 
the center within the fitting uncertainties, only integrated spectra are shown in Fig. HJ We 
give a discussion of the spot morphologies for selected individual sources in Sect. 13.21 

In several sources with complex spectra, there were multiple maser spots for a given 
velocity channel. When these spots are close enough to be blended together, SAD does not 
provide a reliable fit. To overcome this problem, we inspected the fit results for each data 
cube. When multiple spots were blended in the fit results, we repeated the fits manually 
using the task "JMFIT". This procedure is adequate to reliably determine the maser spot 
morphologies for most sources. However, a small number of sources are too complex to 
resolve using this technique. These sources are noted in the discussion in Sect. 13.21 While 
Fig. [T] shows the spot maps determined from our data, higher angular resolution observations 
are required to fully resolve individual spots in these sources and verify our results. 



3.2. Selected individual sources 

G35.03+0.35 - A complex emission structure is seen in this source as expected from its 
line profile. There is an elongated though not perfectly linear structure seen between LSR 
velocities of 41.9 and 44.9 km s~^ with a linear velocity gradient across the structure. In ad- 
dition, there are five other groups of spots covering the f ull velocity range of maser emission. 



The absolute position is consistent with that derived by ICyganowski et al.l (120091 ) although 
their poorer resolution fails to resolve individual spots in channels with multiple spots. The 
24 iim. MIPSGAL image is saturated in the central pixels, and the overall morphology gives 
the appearance of two sources that are blended together. While Table |5] gives the results 
of a single Gaussian fit, higher resolution data is required to determine the presence or ab- 
sence of multiplicity. The source is associated with an extended green object (EGO) in the 
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GLIMPSE data; a more detailed discussion can be found in ICyganowski et al.l ( 120091 ) . 



G35.4O+O.O3 - This is one of three sources for which no rehable counterpart can be de- 
termined at 24 fim. There is a point source 7" away at 18''55™'50'^.87, 2°12'26".4 which 
is connected to extended emission towards the east suggestive of another point source at 
roughly 18'^55'^51''.0, 2°12'20" (no Gaussian fit could be obtained). The latter is about 
3.4" from the maser position. It is not clear whether either of the two sources is associated 
with the 6.7 GHz methanol maser. 



G35.79-0.n- This source was observed by iBartkiewicz et al.l (120091 ) and observed to have 
a linear morphology with a velocity gradient. While we detect this structure in our data, 
there are three other groups of spots, one of which has a ring shap ed morphology. Mo r eover , 



the peak flux density detected in MERLIN is over 20 Jy, while IBartkiewicz et al.l (120091 ) 



detect a peak flux density of only 9.7 Jy. Hence, the EVN observations suffer from missing 
flux, and the overall morphology of the maser emission is complex rathe r than linear. We 
also n ote a significant position offset between the coordi nates reported by IBartkiewicz et al. 



(120091 ) and those in Table [31 However, the position of IBartkiewicz et al. 



mm is derived 



from MERLIN single baseline data which has a much larger uncertainty. Th e posit ion offset 
in this source is consistent with the uncertainties quoted in lSzymczak et al.l (120071 ). 



G 36.70+0.09 - The s pot n iorphology as seen by MERLIN is very similar to that derived 
by IBartkiewicz et al.l (120091 ) using the EVN. There are two main groups of maser spots in 
this source. The southern group shows an inverted "S" structure, while the northern group 
hosts a linear distribution with a velocity gradient between 61.9 and 62.4 km s"^. The EVN 
observations, being more sensitive, also see more maser spots at high velocities (> 63 km s~^), 
while we detect only a single spot at 63.0 km s~^. 

G36.84-0.02 - There are two groups of maser spots in this source (Fig. [2]) that are separated 
by 1.2". However, both groups show spectral features within the overall velocity range 
containing emission. This suggests that the two groups are discrete sources rather than 
being kinematic features of a single source. However, given that the distance to the source 
is only 3.7 kpc, this would imply a physical separation of only 0.02 pc (~ 4500 AU). High 
angular resolution mid-infrared and submillimeter observations are required to determine 
whether the two groups of spots are indeed discrete, or whether they are excited by a single 
source. 

G36. 92+0.48 - The MIPS 24 /zm image shows a cluster of at least 2-3 sources within 10" in 
this region. Taking into account the large distance to this source (15.8 kpc), the multiple 
24 yum point sources probably reside in the same molecular cloud. 



G37.02-0.03 and G37.04-0.04 - G37.04-0.04 is a new individual source detected in the 
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MERLIN data, and is located 49" away from G37.02-0.03. The spectral features of the 
source can be seen at a weak level in the published spectrum of G37.02+0.03 in Paper I, 
and was most probably missed in the original survey due to its close proximity to G37.02- 
0.03. The maser emission in G37.02+0.03 is confin ed to a small region, very similar to 
that observed with EVN by iBartkiewicz et al.l ( l2009l ). We note however that the peak flux 
density as seen in MERLIN is much larger than that in EVN, and that maser spots outside 
the central core appear to be resolved out in the EVN data. In contrast, the morphology of 
G37.04-0.04 is more complex with more spectral features. There are two parallel emission 
features seen, each with a linear morphology. While the northern feature shows a linear 
velocity gradient, the southern feature does not show any coherent velocity structure and 
appears to be a mere superposition of different spectral features. 

G37.47-O.il - This is a complex source with several spectral features. The emission features 
between 60.8 and 63.1 km s"^ form a linear structure with a velocity gradient. However, 
the other spectral features are distributed randomly with a general tendency of features 
at lower v elocity to be located to the south. The overall morphology is similar to that 
derived by ICyganowski et al.l (120091 ) using the EVLA although we do not detect the eastern 
spots seen in their map. Consequently, we do not see the double arc structure surmised by 
Cyganowski et al.l ( 120091 ). It should also be noted that sorne cha nnels showed emission in 
multiple spots which were unresolved by ICyganowski et al.l (|2009[ ). 



G37.53-0.il - There are two groups of maser spots in this source separated by 1.4". At 
a distance of 9.9 kpc, this corresponds to a physical separation of 0.07 pc. High angular 
resolution observations at infrared or submillimeter wavelengths are required to determine 
whether the two groups are excited by a single or by two separate massive young stellar 
objects. The MIPSGAL image of this source is completely saturated and hence the position 
of the 24 /im point source cannot be measured. 

G37.55+0.19 - There are two well separated groups of spectral features in this source which 
are also spatially separated in the spot map (Fig. 3). The features at lower velocity (78.4 to 
80.0 km s~^) are situated about 0.2" to the west of the other spots. The eastern group of 
spots is clustered in accordance with individual spectral features, but show a broad velocity 
trend from the north-west to the south-east with increasing velocity. As in G37.53-0.il, the 
position of the MIPSGAL counterpart cannot be determined due to saturation. 

G37.6O+O.42 - The intensit y in our MERLIN s pectr um is more than a factor of two higher 
than the EVN spectrum of IBartkiewicz et al.l (120091 ). The peak emission in the MERLIN 
data occurs at 87.0 km s~^, while the EVN peak occurs at 85.8 km s~^. Since the masing spot 
at 87.0 km s~^ is about 30 mas to the east and 10 mas to the north of the spot at 85.8 km s~\ 
the reference coordinate quoted in Table [3] is different from that of IBartkiewicz et al.l ( 120091 ). 
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While the spot morphology determined from our data shows similarities to that derived 
from EVN, some differences exist. The MERLIN spots are much more tightly concentrated 
in declination between 84.9 and 89.0 km s~^, forming tw o linear structures with monotonic 
velocity gradients. The two easternmost maser spots of iBartkiewicz et al.l ( 120091 ) at veloc- 
ities lower than 85 km s~^ are not detected in our data. However, we detect several spots 
at velocities greater than 91 km s~^ which are not detected in the EVN data. Another 
qualitative difference is that the spot distri bution between 52 and 56 km s~^ resembles a 
ring rather than the linear structure seen in IBartkiewicz et al.l (120091 ) . 



G37. 74-0.12- The MIPS GAL image of this source shows significant extended emission with 
the background subtracted image revealing a bow-shock morphology. It is not clear whether 
the morphology is intrinsic or whether it is due to a superposition of multiple sources. The 
6.7 GHz methanol maser is coincident with the brightest emission region. The position 
reported in Table [5] is the result of a Gaussian fit to the bright emission region. 

G37. 76-0.19- The morphology of this source as seen with the EVLA shows a compact group 
of spots between 54.5 and 55.1 km s~^ (which is the strongest spectral feature in the source), 
while the other spectral features are distributed over a larger area to the east. The latter 
being weak, have relatively large fitting uncertainties, and so we do not attempt to deduce 
any morphologies from the data. 

G37. 77-0.22- This methanol maser is embedded in a large area of extended 24 //m emission, 
and no point source can be discerned at the maser position. 

G38. 03-0.30 - As with G37.6 the MERLIN spectrum is much stronger than the EVN 
spectrum of IBartkiewicz et al.l ( 120091 ) . and the velocity of peak emission in the MERLIN 
data (58.2 km s^^) is different from that in the EVN data (55.7 km s~^). This leads to a 
small offset between the position quoted in Table [3] and that of IBartkiewicz et al.l ( 120091 ). 
However, the spot morphology determined from MERLIN is consistent with (and essentially 
identical to) that determined using EVN, except for the detection of spots at high LSR 
velocities (62.8 and 63.7 km s"^) in the MERLIN data. 

G38. 12-0.24 ^ There are three groups of spots in this source. The velocity components 
between 68.1 and 71.4 km s~^ are tightly concentrated, while the two other groups of features, 
76.7 - 77.7 km s~^ and 79.1 - 79.3 km s~^, lie to the southeast and northeast respectively. 
The weaker spectral feature around 74 km s~^, detected in only one channel, lies close to 
the central cluster of spots. 

G38.20-0.08 - This is a complex source with several strong spectr al features. We measure a 
peak flux density of 9.7 Jy at 84.3 km s~^, while the EVN data of IBartkiewicz et al.l (|2009[ ) 
detect a peak of only 0.83 Jy at 79.6 km s~^. Paper I lists a single dish peak flux density of 
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11.1 Jy at 79.6 km 3"^ The MERLIN spectrum is very similar to that obtained at Arecibo 
except that the 84 km s~^ feature is much stronger in the MERLIN observation. While this 
can be attributed to variability, the significant discrepancy between EVN and MERLIN is 
unlikely to be due to variability, but suggests missing fiux. The spot distribution is extremely 
complex with some channels showing as many as 3 spots some of which are blended together. 
H ence, the spot map shown in Fig. [1] should be treated with caution. While most of the spots 



of iBartkiewicz et al.l (120091 ) can be identified in our spot map, we also detect several features 



that are not seen in EVN. We also do not detect any maser spot ~ 180 mas north-west of 



the central concentration of spots that is seen in the 



VN data. The absolute position in 



Table [3] is consistent with that of IBartkiewicz et al.l (120091 ) when the different velocity of 
peak emission is taken into account. 

The 24 /im MIPSGAL image shows a strong point source at 19''01"n8".86, 4°39'26".7 
with the methanol maser being coincident with a weak source which lies on the Airy diffrac- 
tion ring of the strong source. The weak source, which cannot be fit by a Gaussian, has an 
approximate position of 19'^01™'18''.7, 4°39'37" by visual inspection. Further work including 
PSF subtraction of the strong source is required to determine the position of the counterpart 
accurately. The position of the MIPS counterpart reported in Tabled] is obtained by visual 
inspection. 

G38.26-0.08 - This source (observed with EVLA) has three groups of spectral features 
which are spatially resolved. The features 6.3 - 7.6 km s~^ and 11.7 - 12.9 km s~^ populate 
the eastern and south-eastern part of the spot map respectively, while the strong spectral 
features between 14.2 and 15.5 km s~^ are more compact and lie to the west. 

G4O.28-O.22 - This is another complex source with an unusually large number of spectral 
features. The frequency used to observe this source led to the loss of features at velocities 
greater than ~ 79 km s~^. The overall shape of the spectrum obtained using MERLIN is 
very similar to that in Paper I although variability is seen in some spectral features - the 
feature around 65 km s~^ is stronger while the peak fiux density is lower in the MERLIN 
data compared to Paper I. The spot map shows about 8 groups of spots, at least two of which 
are extended linearly with a monotonic velocity gradient. However, the overall distribution 
of spots is complex with no velocity structure. The 24 /xm MIPSGAL image is saturated in 
the center. 

G4O.62-O.I4 ^ There are four spectral features in this source that are clearly resolved spa- 
tially. However, the overall morphology should be classified as "multiple" if not complex. 
The MIPSGAL 24 //m counterpart is fully saturated and hence its position could not be 
measured. 
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G41-12-0.22- There are two spectral features in this source with the stronger feature showing 
a hnear morphology with a velocity gradient. 

G41-23-0.20 - There are four groups of spots in this source. While three groups show 
linear morphologies, two out of three groups do not show any clear velocity gradients. The 
MIPSGAL 24 /im image shows extended emission at the maser location. The background 
subtracted image shows a few point sources along with some residual extended emission. 
Hence there is some uncertainty in the position of the counterpart reported in Table O 

G41-34-0.14 ^ The maser spots in this source fall into two broad groups. The spots at 
velocities less than 10 km s~^ lie to the southeast, while those at velocities greater than 
10 km s~^ lie to the northwest. No distinct morphology can be discerned in the northwest, 
while the southeast group could possibly be fit with a ring morphology. 

G42.O3+O.I9 - This is a complex source. The emission between 7.2 and 11.0 km s~^ are 
distributed in four clusters with an overall north-south trend with decreasing velocity. The 
other spots are distributed in a complex morphology though it is conceivable to fit the 
spots between 7.2 and 14.6 km s~^ with a ring distribution. We also detect a maser spot at 
17.0 km s~^ although the MERLIN spectrum does not show a spectral feature at this velocity 
(this spectral feature is clearly seen in the Arecibo spectrum in Paper I). The position of 
the 24 ixm. counterpart has some uncertainty since a significant fraction of the PSF core is 
saturated and blanked. 

G42.43-O.26 - The position of the 24 /im counterpart of this maser cannot be measured on 
account of it being completely saturated. 

G42.7O-O.i5 - There are five spectral features in this source each of which are spatially 
distinct. While most individual features display a linear or arched morphology, there is no 
overall velocity structure in this source. It is possible that the different maser spots can be 
fit with a ring morphology. 

G43.O4-O.46 - There are three regions of emission in this source, two in the southeast and 
one in the northwest. The northwest group is offset from the other two groups by ~ 2.5", 
which corresponds to a physical separation of 0.1 pc. While it is possible that the northwest 
group of spots is a distinct maser excited by a different YSO than that pumping the southeast 
group, complementary data at high angular resolution is required to confirm this scenario. 
The MIPSGAL counterpart is close to the southeast group of spots. 

W49N region - Using the MERLIN data, we are able to obtain accurate astrometry for all 
five methanol masers identified in Paper I. We do not identify any new emission regions to 
a la limit of 60 mJy beam~^. Four out of five sources have previously published positions 
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accurate to 0.4" using the Australia Telescope Compact Array (ATCA; ICaswellll2009l ). The 
two weak sources, G43. 17-0.00 and G43. 18-0.01 show a simple compact morphology or are 
detected in only one maser spot. The other three sources are discussed in more detail 
below. The entire region is saturated in the MIPSGAL 24 /xm image and hence no 24 /im 
counterparts can be determined for any of these sources. 

G43. 15+0. 02- Although the spectrum in Fig. [T]is affected by negative sidelobes of G43. 16+0.02, 
the spot map was verified to not suffer from these effects. The spot map shows a simple 
linear morphology with an overall velocity gradient. 

G43. 16+0.02 - The spectral features between 7.7 and 8.9 km s~^, and 9.1 and 10.2 km s~^ 
form two roughly linear features with monotonic velocity gradients (though a sub-feature 
around 8.9 km s~^ is an outlier) that are perpendicular to each other. The emission between 
15 and 20 km s~^ traces a linear and ring shaped morphology in addition to a few relatively 
isolated spots. 

G43. 17+0.01 - There are at least five morphological features in this source. The emission 
between 18.8 and 19.3 km s~^, and 20.0 and 20.5 km s~^ traces two linear structures in the 
northwest, while the 20.2 - 20.7 km s~^ and 21.2 - 22.1 km s~^ features form two clusters 
to the southeast. The emission between 20.2 and 21.1 km s~^ traces another linear feature 
which connects to the cluster between 21.2 and 22.1 km s~^. In addition, a number of spots 
between 19.1 and 20.0 km s~^ form a relatively random distribution. 

G43.8O-O.I3 - This is a complex source with several maser spots being blended in our 
MERLIN data. While the spot distribution in Fig. [1] shows two linear features between 39 
and 41 km s~^, higher angular resolution data is required to confirm this. We also detect 
four other clusters of spots in this source. The spectrum in Fig. [T] is Banning smoothed to 
suppress ringing arising from the Gibbs phenomenon. The peak flux density is 53.9 Jy in the 
unsmoothed spectrum which reduces to 40.8 Jy after Hanning smoothing. Considering the 
velocity resolution of the MERLIN data, this source is consistent with no variability seen 
between the epochs of the Arecibo and MERLIN observations. The position of the maser's 
24 fim counterpart cannot be measured on account of saturation. 

G44-31+0.04 - The 24 //m counterpart of this maser is saturated in the central pixels. 

G45.O7+O.I3 - Even though the peak flux density is over 60 Jy in this source, there are 
only two spectral features with the strong feature displaying a linear morphology with a 
monotonic velocity gradient. The MIPSGAL 24 fim image is fully saturated at the location 
of the maser, and so no point source counterpart could be identified. 

G45.44~l'0.07 - This source is in a bright region with extended emission at 24 jj,m that is 
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saturated close to the maser position. Although no point source is seen after background 
subtraction, this could be due to the close proximity to blanked region. 

G45.47+O.I3 - There are four emission regions in this source. Of particular interest is the 
double emission structure at 59.6 and 59.7 km s"^ (i.e. each velocity channel has emission 
in two spots) in the northeast. As with G45.44+0.07, this maser is located at the edge of a 
saturated region at 24 /im, and hence no 24 /im can be identified for the source. 

G45.47+O.O5 - While there are multiple spectral features in this source, the maser spots 
between 55.6 and 57.6 km form a single morphological feature although the velocity 
gradient is not monotonic across the spots between 56.5 and 57.6 km s~^. A second linear 
feature between 57.7 and 58.3 km s~^ is located about 10 mas to the east. 

G45.49+O.I3 - The maser spots in this source are concentrated in a compact region. The 
spot morphology is linear with a velocity gradient, except for the spot at 58.0 km s~^. 

G45.8I-O.36 - There are five emission regions in this source. The spectral features between 
velocities of 59.7 and 66.4 km s~^, which include the strongest peak, each display a simple 
morphology. In contrast, the features between 67.4 and 70.3 km s~^ show a more complex 
distribution with no overall velocity structure. 

G46. 12+0.38 - This source displays three spectral features, two of which have linear mor- 
phologies. The third feature is identified in only two spots, and hence no morphological 
information can be deduced for this feature. 

G48.9O-O.27 - While the overall morphology in this source is linear, there is no velocity 
structure across the feature. 

G48.99-0.30- This maser, as with a number of other sources in the W51 region is in a region 
that is saturated at 24 //m and hence no point source counterpart can be identified at this 
wavelength. 

G49.27+O.3I - There are two emission regions in this source, one to the northeast and the 
other to the southwest. The southwest region shows a single extended morphology between 
velocities of -6.1 and -1.9 km s~^, and a monotonic velocity gradient is seen between -5.6 
and -3.0 km s~^. 

G49.35+O.4I - The emission in this maser is relatively compact with the spots between 
velocities of 66.9 and 69.1 km s~^ forming a reversed "S" shaped morphology. 

G49.4I+O.33 & G49.42+O.32 - The single dish spectrum shows two main groups of spec- 
tral features, one between -27 and -23.5 km s^^, and a second group between -16 and 
-9.5 km s~^, although Paper I also shows weak emission around -19.5 km s~^. As indicated 



- 16 - 



by ICyganowski et al.l (|2009[ ). at high angular resolution, the two groups of features separate 
into two distinct regions that are separated by ~ 8" (0.47 pc physical separation at a distance 
of 12.2 kpc). The two sources are hence treated as distinct masers. Within our la limit of 
45 mJy beam"^, the western source (049.41+0.33) only shows emission between -16 and 
-10 km s~^, while the eastern source (G49.42+0.32) only shows emission between -27 and 
-23.5 km s~^. We thus do not detect any maser s pots a t velocities less than -23 km s~^ 
in G49.41+0.33 as was seen by ICyganowski et al.l (120091 ). The high angular resolution of 
MERLIN clearly resolves the morphologies of both sources. G49.41+0.33 consists of a set of 
four morphological features extending in the east-west direction. Each feature is extended 
in the north-south direction as shown in Fig. [H In contrast, G49.42+0.32 consists of five 
groups of maser spots spanning 140 mas in the north-south direction. 

While G49.42+0.32 has a clear 24 fim counterpart, the case of G49.41+0.33 is am- 
biguous. There appears to be multiple point sources in the Airy ring of the source associ- 
ated with G49.42-f0.32, one of which is about 3" away from G49.41+0.33 (coordinates of 
19''20™59^0, 14°46'52" by visual inspection). Although it is not clear whether G49 .41+0.33 
has a 24 /im counterpart, it is associated with an EGO (lOyganowski et al.ll2009l ). Higher 
spatial resolution mid-infrared observations are required to determine the counterpart of 
G49.41+0.33. 

W51 main - There are four methanol masers in the W51 main region identified in Paper I. 
This region was n ot observed in t he present observing program. However, we have used the 
MERLIN data of IXu et al.l (120091 ) to determine the positions of each of the four sources in 
the channels of strongest emission (see Table [3]). The dynamic ra nge of the data is limited 
by the very strong source G49.49-0.39 (W51e2 in IXu et al.l 120091 ). Moreover, the velocity 
resolution of this data is only 0.72 km s~^ after Hanning smoothing (required to reduce the 
Gibbs ringing from G49.49-0.39). Hence, no spot maps are d erived fo r any of the sources. 
Our pos i tions are in good agreement with those determined by ICaswelll (120091 ) using ATCA. 
Caswelll (|2009l ) also identified a n additional site of emission in thi s region. Several regions of 
emission were also detected by iPhillips fc van Langeveldd (|2005l ) using the EVN, although 
they do not provide any absolute positions for the sources. The entire region is saturated at 
24 /im and hence no 24 /im counterparts can be determined for any of the sources. 

G49.6O-O.25 - This is an extremely complex source with several emission channels having 
multiple maser spots that are blended together. Although the spot map in Fig. [1] is produced 
by manually fitting synthesized beams to the blended spots, observations at higher angular 
resolution are required to verify the same. The emission morphology is also extremely 
complex. The MIPSGAL 24 /im image suggests the presence of two sources that are blended 
together, with another source located on the Airy ring. The counterpart position indicated 
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in Table |5] was obtained by fitting two Gaussians to the central source(s). 

G52. 92+0.41 - This is an interesting source with three spectral features forming a linear 
north-south morphology. In addition, each spectral feature also shows a north-south mor- 
phology with a velocity gradient. However, there is no overall velocity gradient across the 
three features. The 24 fim counterpart is saturated at the central pixels. 

G53.O4+O.II - The 24 yum image shows a cluster of at least four sources with the central two 
sources blended with each other. The counterpart listed in Table E] was obtained by fitting 
two Gaussians to the central pair of sources. 

G53. 14+0.07 - No 24 nm counterpart could be identified for this maser due to it being 
saturated. 

053. 62+0. 04 - The central pixels of the maser counterpart in the 24 /xm MIPSGAL image 
are saturated. 



Discussion 



4.1. Comparison with previous work 

Of the 82 sources with good astrometry, six (G35.79-0.17, G36.80+0.09, G37.02-0.03, 
G37.4 7-0.il, G38.03-0.30 and G38.20-0.08) have been studied using the EVN bv lSartkiewicz et al, 
J2009h . five (0 35.03+0.35, G37.4 7 -0.il , G37.53-0.il, G49.41+0.33 and 049.42+032) using 



the EVLA by ICyganows 
(see Table 1 of Oaswell 



d et all mm . and fourteen using the ATCA by ICaswelll (120091 ) 
20091 ). In addition, the absolute positions of 043. 15+0. 02 and 



O49.41+0.33 (IRAS 19186+ 1440), der i ved fr om MERLIN (but at poorer velocity resolu- 
tion) have been published by IXu et al.l (120091 ). For the EVN sources, most of our positions 
agree with previous work within our MERLIN astrometric uncertainty of 15 mas. The ex- 
ceptions are 035. 79-0. 17 and G37.47-0.il which have position differences of 140 mas and 
270 mas respectively. However, the positions for these two sources wer e determined from 



MERLIN single baseline data (no EVN positions derived; see Sect. 3.1 of iBartkiewicz et al. 



20091 ) and have unc ertainties of 175 raas and 365 mas respectively (ISzymczak et al.l 120071 ). 
The EVLA data of ICyganowski et al.l (120091 ) were taken in the B configuration, for which 
uncertainties larger than 100 mas can be expected for the absolute astrometry. Taking this 
into consideration , there is good agreement between the positions in Table |3] and those of 



Cyganowski et al. 



(I2OO9I). The positions derived from ATCA have an root mean sq uare un- 



certai nty of 0.4" (ICaswellll2009l ). Our positions in Table [3] agree well with those of lOaswell 



(120091 ) except for 043. 17-0.00 and 43.80-0.13. The discrepancy in the case of G43.17-0.00 is 
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due to a typing mistake in the right ascension of ICaswelll ( 120091 ) (the right ascension should 
be 19^10'"16'^.72 instead of 19^10''"16'^.2 7; Caswell, priv ate communication). There is also 
excellent agreement with the positions of lXu et al.l ( 120091 ). with minor differences presumably 
arising from the poorer velocity resolution of their data (which can lead to blended spots). 

Comparing the positions in Tables [3] and H] with those of Paper I, we find a mean offset 
between the accurate position and the single dish position (using a Gaussian probability 
density function) to be -8.5" (-0^.56) in right ascension and 1.2" in declinatio n. This confirni s 
the observation of a bias in the right ascension of the Arecibo positions by ICaswelll ( l2009l ). 
However, the cause of this offset is not clear. It is unlikely to be due to telescope pointing 
since no bias is seen in the pointing scans in the original data. 



4.2. Spatial extent of maser emission 

When analyzing the spatial extent of maser emission, we exclude sources that have sus- 
pected multiple YSOs exciting well separated spatial groups of spots (G36.84-0.02, G37.53- 
0.11 and G43.04-0.46) in addition to sources that have fewer than four maser spots. We 
derived spatial extents using the following procedure. We first defined the extent of the 
maser spot distribution in right ascension and declination, and carried out a coordinate ro- 
tation to determine an axis along with the extent was minimum. This was defined to be the 
minor axis, with the perpendicular axis being taken as the major axis. The spatial extent 
along the major axis ranges from 26 AU to 7680 AU with a median value of 775 AU, while 
the minor axis ranges from 2 AU to 2400 AU with a median value of 180 AU. Fig. [5] shows 
the distribution of the spatial extent along the major axis. It can be seen that the majority 
of the sources have spot distributions with major axis less than 1000 AU. 

It is of interest to look for correlations between the spatial major axis extent of maser 
emission (the spatial extent along the minor axis is affected by the three-dimensional geom- 
etry of the source) and other parameters such as the velocity width of emission and maser 
luminosity. Fig. IH] shows the spatial extent as a function of the range in velocity of maser 
emission (i.e. AV^ = Vmax — Vmin where Vmin and Vmax are the minimum and maximum 
velocities of maser emission respectively) along with a linear fit to the data. While there 
is significant scatter, it can be seen that sources that have a larger velocity width tend to 
have spots distributed over a larger spatial extent. Fig. [H] clearly shows that the minimum 
value of the spatial major axis extent of maser spots (for a given maser velocity width) has 
a strong dependence on the velocity width of maser emission. 

Fig. [3 shows the spatial major axis extent of maser emission as a function of maser 
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luminosity (taken from Paper III). A power law fit to the data (dashed line in Fig. [7]) gives 
an exponent of 0.51. However, as in Fig. [6l there is significant scatter in this relationship. 
It would be of particular interest to look for correlations between the spatial extent of maser 
emission and the luminosity of the central YSO that excites the maser. However, this 
requires determination of spectral energy distributions of the central object from millimeter 
to mid-infrared wavelengths, which we reserve for a future paper. 



4.3. Methanol maser morphologies 



Figs. [T] and [2] make it clear that it is difficult to classify methanol maser morphologies. 
However, there are a few broad classes: (i) Linear/arched where the overall morphology 
forms a rough straight line or a smooth curve although there may not be any systematic 
velocity gradient; (ii) Compact/simple morphologies where the extent of the maser spots is 
less than ~ 20 mas; (iii) Multiple, where different spectral features are well separated into 
compact clusters of maser spots (the pair morphology used in previous studies would be a 
special case of this class); and (iv) Complex where no regular structure is seen. 

Using the classification scheme above, there are 9 sources with linear /arched morphol- 
ogy, 5 with compact /simple morphology, 12 with multiple groups and 24 with complex 
morphology (see Tables [3] and H]) . It is remarkable that we do not see any source with a clear 
ring morphology although there are a few sources (G35.79-0.17, G36.70+0.09, G41.34-0.14, 
G42.03-F0.19, G42.70-0.15 and G43. 16+0.02) where there are potential ring structures. How- 
ever, these structur es are usually embedded in a larger more complex distribution of maser 
spots. In contrast, iBartkiewicz et al.l (120091 ) found 29% of their sample to have ring like 
emission morphologies. The latter may be a selection effect due to their very high spatial 
resolution resolving a significant fraction of the overall emission. For the five sources in com- 
mon for which emission morphologies have been determined, four have miss ing fiux between 



70 an d 90% . Furthermore, a compar ison between the s i ngle d ish fiuxes of ISzymczak et al. 



(I2OOOI . I2OO2I ) and the EVN fiuxes of IBartkiewicz et al.l (120091 ) shows a similar trend for a 
large number of sources. 

The large fraction of missing fiux in some sources indicates the presence of emission 
structures on scales that are resolved by the E VN. This is rnost lik ely to arise from the 
structure of the maser spot itself. For example, iMinier et al.l (120021 ) found a majority of 
masing spots to consist of a compact core with a halo of extended emission that is resolved 
even on short VLBI baselines. However, if some spots do not have a compact core, they 
would not be detected in VLBI obse rvations. The con i pariso n between the overall emis- 
sion morphology in Fig. [T]and that of IBartkiewicz et al.l (120091 ) suggests examples for both 
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scenarios. For example, the spot morphology in G37.02-0.03 obtained using MERLIN is 
remarkably similar to that of EVN in spite of the latter missing 90% of the flux, suggesting 
that most of the maser spots in this source have a core/halo structure. In contrast, G35.79- 
0.17 and G38.20-0.08 show several maser spots in the MERLIN data that are not seen using 
EVN suggesting that a number of spots in these sources do not have a compact core. Al- 
though a detailed discussion about the structure of a maser spot is beyond the scope of this 
pa per and would require c o mbini ng this data with higher resolution observations, the work 



of iHarvey-Smith fc CohenI (120061 ) shows that in very nearby regions, even MERLIN obser- 
vations could resolve some extended emission. We discuss this in the context of variability 
in Sect. KM 

In light of thi s discussion, it app e ars lik ely that at least some of the ring shaped emission 



structures seen by iBartkiewicz et al.l (120091 ) are embedded in or are part of larger structures. 
This casts doubt on simple kinematic models for such sources, and highlights the dangers 
of interpreting kinematics solely based on morphology. To obtain a better understanding of 
the origin of methanol maser emission in the context of a massive young stellar object, one 
requires either proper motion studies through multi-epoch VLBI, or high angular resolution 
data in mid-infrared and submillimeter wavelengths. 



4.4. Variability of AMGPS methanol masers 

Methanol masers are known to be vari able, and a number of studies have systematically 



studied the variability of the masers (e.g. iGoedhart et al.ll2004j ). To look at the variability 



of the AMGPS masers, we compared the integrated spectra from MERLIN and EVLA with 
those in Paper I. Among the 77 sources with spectra in both MERLIN/EVLA and Paper I, 
23 have similar peak flux densities, 16 have become stronger and 38 are weaker. Among the 
variable sources, 20 sources have changes in their peak flux densities by more than 50%, the 
most dramatic being G35. 59-1-0. 06 which has dropped from 0.82 Jy in the Arecibo data to 
0.20 Jy in the EVLA data. It is interesting to note that all the masers in the W49N region 
have become weaker with three sources being weaker by more than 50%. 

While some MERLIN so urces could suffer from resolving out diffuse emission (e.g. 



Harvey-Smith fc Cohen! l2006l ). this is unlikely for the more distant sources. In any case, 
the poorer resolution of the EVLA is expected to recover all the flux in the source. Among 
the 25 EVLA sources, 5 have the same peak flux densities as Paper I, while 15 are weaker 
(although the change is small in three sources) and 5 are stronger. We however caution that 
this statistic has a bias towards sources that have become weaker since they are more likely 
to be not detected in the MERLIN data. Nonetheless, it appears that a larger fraction of 
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the AMGPS masers have become weaker compared to the flux measurement in Paper I. It 
is interesting to note a similar trend in Paper I when compared with previous work (see 
Table 2 of Paper I). 



4.5. Association with mid-infrared emission 



The association between mid-infrared emission and 6.7 GHz methanol masers has been 
explored in a number of previous works. While early studies found a nuniber o f sources to not 
have mid-infrared counterparts at 10.5 /im and/or 20.0 /xm ( iWalsh et al.ll200ll ) or in the MSX 
point source catalog (e.g. Paper II), recent work has found methanol masers to be generally 
associated with mid-infrared emission. For example, lEllingsenI ( l2006l ) found most methanol 
masers to have a count e rpart in the Spitze r GLI MPS E survey within 5^^ . The more recent 
work of IXu et al.l (120091 ) . ICyganowski et al.l (120091 ). and iBartkiewicz et al.l (120091 ) have found 
very close correspondence between 6.7 GHz methanol masers and mid-infrared emission at 
24 fim and 4.5 fim. 

Fig.|8]shows the histogram of the angular physical separation between 6.7 GHz methanol 
masers and their 24 ^um counterparts for the 60 sources in Table [51 The median angular sep- 
aration between a methanol maser and its 24 fim counterpart is 0.9" and the corresponding 
physical separation is 7100 AU. We also note that 18 sources had no counterpart identi- 
fied due to detector saturation, and only four sources had no counterpart identified in the 
MIPSGAL data. Of the four, two were in regions of extended emission, and an additional 
source has a potential counterpart on the Airy ring of a stronger source. Hence, almost all 
6.7 GHz methanol masers have associated mid- infrared emission. These results also show 
that the lack of counterparts in the MSX Galactic plane survey is primarily a result of poorer 
sensitivity and to a lesser extent the poorer resolution of the MSX satellite. The 24 /xm point 
so urce fluxes of the me thanol maser counterparts are frequently less than 1 Jy (see Table 3 
of jPandian et al.l |2010| ) which would be undetected by MSX at 21 /im, while the poorer 
resolution of MSX would lead to significant confusion in some crowded fields. 

This close association of mid- infrared point sources associated with 6.7 GHz methanol 
masers is important in the context of theoretical models for the maser emission. For a wave- 
length of 24 /zm, the Planck function has its maximum at 120 K, which is very close to the 
temperature of the molecular hot cores that the high mass protostellar objects exciting class 
II methanol masers reside in. This wavelength is also close to that of one the main bands 
invoked in radiative class II meth anol maser pumping schern es (excitation from the torsional 
ground state, Vt = to Vt = 2; see lSobolev fc Deguchilll994l ). These two facts together, com- 
bined with the ubiquitous presence of 24 /im sources apparently ensure extremely effective 



- 22 - 



excitation of the masers. 



The work of iBartkiewicz et aLl (|2009[ ) found that masers with ring shaped and regular 
morphologies had strong 8 ^m emission. Particularly, they found that no source that was 
saturated at 8 /xm had a complex morphology. While we have not inspected the 8 fim 
images of the AMGPS sources, we do not see a similar clear correlation with the 24 /im data. 
Considering the 16 sources that were saturated at 24 /im (even if only in the central pixels) 
and had morphologies determined, there were 5 linear /arched, 3 multiple, one compact and 

7 complex sources. There is no statistically significant difference between the morphologies 
of masers associated with sources that are saturated at 24 /xm and that of the overall sample. 
However, further work is required to verify whether a similar trend is seen with the GLIMPSE 

8 //m data. 

It is interesting to note that from Fig. [HI the separation between the maser and mid- 
infrared source does not peak at zero. This could arise either from the uncertainty in the 
absolute astrometry of MIPSGAL, or from astrophysical reasons u nderlying the location of 



the maser emission relative to the source of mid-infrared emission. ICarey et al.l ( l2009l ) ana- 
lyze the relative positions of 24 fim MIPSGAL point sources with those of 8 /im counterparts 
in the GLIMPSE point source catalog, and find the median astrometric offset between the 
MIPSGAL and GLIMPSE to be 0.85". This is very close to the median separation between 
methanol masers and MIPSGAL counterparts, and may be the likely explanation for this 
phenomenon. We have not undertaken the exercise of cross-correlating the maser positions 
with the GLIMPSE catalog due to the presence of emission and absorption features from dust 
grains in so me IRAC bands, wh i ch oft en leads to extended emission. Moreover, it is clear from 



the work of ICyganowski et al.l ( l2009l ) that many infrared counterparts of methanol masers 
are also extended at 4.5 /xm. We thus reserve addressing the relation between GLIMPSE 
point sources and the AMGPS masers for a future paper. However, the issue of whether 
or not there is a systematic offset between mid-infrared sources and methanol masers for 
astrophysical reasons can be resolved by high spatial resolution imaging from ground based 
facilities. 

An interesting consequence of the close correspondence between MIPSGAL 24 fim point 
sources and methanol masers is that one can use the former to estimate positions for masers 
that have not been detected with an interferometer. In the AMGPS catalog, there are six 
sources which are not detected, or haven't been observed due to constraints of telescope time. 
Since the la uncertainty in the AMGPS position is 7", we looked for 24 //m point sources 
within 18" of each of the six sources. For five sources, there is a clear 24 /jm counterpart, 
while for G36.90-0.41, there are two MIPSGAL point sources within 18" of the Arecibo 
position. The coordinates of the 24 /xm sources are listed in Table [61 Considering the 
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statistics of the offset between MIPSGAL point sources and metfianol masers, we suggest 
that these positions can be used to refine the single dish positions to an accuracy of ~ 1". 

5. Conclusions 

We have determined the absolute positions of 82 methanol masers in the region surveyed 
by the AMGPS including two sources that were not present in the original catalog. The me- 
dian spatial extent along the major axis of maser emission is ~ 775 AU. The spatial extent 
is correlated to both the velocity width of the maser and the maser luminosity though there 
is significant scatter in both relationships. Among the 50 sources for which the morphologies 
of maser emission could be determined, 9 sources have a linear/arched morphology, 5 have 
compact/simple structures, 12 have multiple compact groups of emission and 24 have com- 
plex morphology. In sources where there is a large discrepancy between the flux recovered in 
VLBI experiments and that measured with a single dish, we often find the missing fiux to be 
distributed over complex structures. Hence, some relatively simple morphological structures 
seen in VLBI are likely to be part of a more complex emission structure in at least some 
sources. This makes interpretation of the origin of maser emission difficult, and multi-epoch 
proper motion studies or complementary high resolution mid-infrared or submillimeter obser- 
vations are required to resolve this problem. We also find very close correspondence between 
mid-infrared emission and 6.7 GHz methanol masers as would be expected from theoretical 
models of the maser emission. 

We are extremely grateful to Anita Richards for invaluable help in MERLIN data pro- 
cessing and other useful discussions. We also thank Rob Bcswick and Peter Thomasson for 
their assistance in MERLIN data reduction, Anna Bartkiewicz for useful discussion, and 
Jim Caswell whose comments helped improve the paper. YX was supported by the Chinese 
NSF through grants NSF 11073054, NSF10673024, NSF 10733030, NSF 10703010 and NSF 
10621303, and NBRPC (973 Program) under grant 2007CB815403. This work has benefited 
from research funding from the European Community's sixth Framework Programme under 
RadioNet R113CT 2003 5058187. This work was carried out in part by the Jet Propul- 
sion Laboratory, California Institute of Technology. This research has made use of NASA's 
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Table 1. Positions of the phase cahbrators used in the MERLIN observations. 



Source name 


RA (J2000) 


Dec. (J2000) 






Reference 




(h m s) 


(0 , 


(mas) 


(mas) 




B1904+013 


19 07 11.996253 


01 27 08.96250 


1.12 


1.35 


1 


B1854+061 


18 56 31.838813 


06 10 16.76649 


0.62 


0.87 


2 


B1919+086 


19 22 18.633756 


08 41 57.37052 


1.80 


4.03 


3 


B1920+154 


19 22 34.699314 


15 30 10.03262 


0.73 


0.78 


4 


B1922+155 


19 24 39.455870 


15 40 43.94172 


0.32 


0.26 


4 



Note. — The columns show the source name, equatorial J2000 coordinates, un- 
certainties in right ascension and dechnation, and the reference for the same. 

References. — (1) Petrov et al. 2005; (2) Fomalont et al. 2003; (3) VLBA 



cahbrator manual: http://www.vlba.nrao.edu/astro/calib/vlbaCalib.txtl (4) Beasley 
et al. 2002 



Table 2. Positions of the phase calibrators used in the EVLA observations. 



Source name 


RA (J2000) 


Dec. (J2000) 


(Tq 




Reference 




(h m s) 




(mas) 


(mas) 




J1851+005 


18 51 46.723009 


00 35 32.36604 


1.82 


2.30 


1 


J1856+061 


18 56 31.838813 


06 10 16.76649 


0.62 


0.87 


1 


J1922+155 


19 22 34.699314 


15 30 10.03262 


0.73 


0.78 


2 



Note. — The columns show the source name, equatorial J2000 coordinates, un- 
certainties in right ascension and declination, and the reference for the same. 



References. — (1) Fomalont et al. 2003; (2) Beasley et al. 2002 



Table 3. Absolute positions of the AM GPS 6.7 GHz methanol masers determined using 

MERLIN. 



Source name RA (J2000) Dec. (J2000) I b Sp Vp V^ys Morphology 





(h m s 


) 


(°' 


") 






(°) 


(°) 


(Jy) 


(km s-i) 


(km s-i) 




G35.03+0.35 


18 


54 


00, 


.658 


02 


01 


19 


.23 


35, 


.02469 


0.34981 


17.8 


44.4 


52.8 


C 


G35.25-0.24 


18 


56 


30, 


.388 


01 


57 


08 


.88 


35, 


.24726 


-0.23675 


0.90 


72.4 


62.0 


L/A 


G35.79-0.17 


18 


57 


16, 


.892 


02 


27 


58, 


.05 


35, 


,79267 


0.17456 


22.2 


60.7 


60.0 


M 


G36.64-0.21 


18 


58 


55, 


.236 


03 


12 


04, 


.72 


36, 


,63362 


0.20287 


0.74 


77.5 


75.2 




G36.70+0.09 


18 


57 


59, 


.123 


03 


24 


06, 


.12 


36, 


.70526 


0.09630 


7.5 


53.1 


59.7 


C 


G36.92+0.48 


18 


56 


59, 


.786 


03 


46 


03 


.60 


36, 


.91820 


0.48293 


0.47 


-35.9 


-30.7 




G37.02-0.03 


18 


59 


03, 


.642 


03 


37 


45, 


.08 


37, 


.03020 


-0.03847 


6.9 


78.4 


80.6 


S 


G37.04-0.04 


18 


59 


04, 


.406 


03 


38 


32 


.77 


37, 


,04343 


-0.03524 


6.0 


84.7 


80.6 


c 


G37.47-0.il 


19 


00 


07, 


.144 


03 


59 


53 


.08 


37, 


.47897 


-0.10467 


5.8 


58.3 


58.8 


c 


G37.53-0.il 


19 


00 


16, 


.056 


04 


03 


16, 


.09 


37, 


.54606 


-0.11182 


2.5 


50.1 


53.0 


M 


G37.60+0.42 


18 


58 


26, 


.799 


04 


20 


45, 


.47 


37, 


,59784 


0.42522 


16.3 


87.0 


89.5 


c 


G38.03-0.30 


19 


01 


50, 


.470 


04 


24 


18, 


.94 


38, 


.03748 


-0.30019 


8.3 


58.2 


62.0 


c 


G38.12-0.24 


19 


01 


44, 


,152 


04 


30 


37 


.42 


38, 


,11894 


-0.22873 


3.0 


70.7 


83.2 


M 


G38.20-0.08 


19 


01 


18, 


.730 


04 


39 


34, 


.32 


38, 


,20322 


0.06660 


9.7 


84.3 


83.3 


c 


G38.66+0.08 


19 


01 


35, 


.244 


05 


07 


47, 


.36 


38, 


,65275 


0.08767 


1.3 


-31.3 


-39.2 




G38.92-0.36 


19 


03 


38, 


.659 


05 


09 


42, 


.49 


38, 


.91576 


-0.35296 


2.4 


32.1 


38.7 


s 


G39.39-0.14 


19 


03 


45, 


.312 


05 


40 


42, 


.68 


39, 


.38764 


-0.14060 


0.26 


60.4 


66.1 




G40.28-0.22 


19 


05 


41, 


.215 


06 


26 


12 


.69 


40, 


,28179 


-0.21929 


20.5 


74.1 


73.0 


c 


G40.62-0.14 


19 


06 


01, 


.630 


06 


46 


36 


.18 


40, 


.62249 


-0.13821 


17.4 


31.1 


32.3 


M 


G40.94-0.04 


19 


06 


15, 


.378 


07 


05 


54, 


.49 


40, 


.93438 


-0.04086 


3.4 


36.6 


40.2 


M 


G41. 12-0.22 


19 


07 


14, 


.856 


07 


11 


00, 


.69 


41, 


.12299 


-0.22011 


1.3 


63.6 


60.1 


M 


G41.23-0.20 


19 


07 


21, 


.378 


07 


17 


08, 


.17 


41, 


,22603 


0.19708 


5.1 


57.2 


59.2 


M 


G41.34-0.14 


19 


07 


21, 


,842 


07 


25 


17 


.27 


41, 


.34753 


-0.13626 


11.8 


11.7 


12.6 


C 


G42.03+0.19 


19 


07 


28, 


.185 


08 


10 


53 


.47 


42, 


.03438 


0.19025 


16.1 


12.8 


18.0 


c 


G42.30 0.30 


19 


09 


43, 


,592 


08 


11 


41, 


,41 


42, 


,30354 


0.29923 


6.2 


28.1 


27.1 


c 


G42.43-0.26 


19 


09 


49, 


,858 


08 


19 


45, 


.40 


42, 


,43474 


0.26011 


1.6 


66.8 


64.5 




G42.70-0.15 


19 


09 


55, 


.069 


08 


36 


53 


.45 


42, 


.69798 


-0.14735 


3.2 


-42.9 


-44.3 


c 


G43.04-0.46 


19 


11 


38, 


.984 


08 


46 


30 


.71 


43, 


.03797 


-0.45282 


8.7 


54.8 


57.3 


M 


G43.08-0.08 


19 


10 


22, 


.050 


08 


58 


51, 


,49 


43, 


,07404 


-0.07686 


6.5 


10.2 


12.3 


C 


G43.15+0.02 


19 


10 


11, 


,049 


09 


05 


20, 


,46 


43, 


,14897 


0.01319 


10.2^ 


13.3 


10.9 


L/A 


G43.16+0.02 


19 


10 


12, 


.883 


09 


06 


12 


.22 


43, 


.16521 


0.01314 


20.1 


9.3 


3.5 


C 


G43.17+0.01 


19 


10 


15, 


.356 


09 


06 


15 


.50 


43, 


.17072 


0.00453 


5.9 


20.2 


12.0 


c 


G43.17 0.00 


19 


10 


16, 


.720 


09 


05 


51, 


.27 


43, 


,16734 


-0.00355 


1.6 


1.1 


1.7 


s 


G43.18-0.01 


19 


10 


20, 


.068 


09 


05 


55 


.87 


43, 


.17483 


-0.01518 


0.32 


11.2 


11.9 




G43.80-0.13 


19 


11 


53, 


.990 


09 


35 


50, 


.61 


43, 


.79542 


-0.12709 


53.9'' 


39.6 


43.7 


c 


G45.07+0.13 


19 


13 


22, 


.129 


10 


50 


53 


.11 


45, 


.07122 


0.13213 


63.5 


57.8 


59.2 


L/A 
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Table 3 — Continued 



Source name 


RA (J2000) 


Dec 


. (J20G0) 




; 

L 





c 

Op 


T/ 


T/ 

^sys 


JVlorpnology 




(h m s) 


(°' 


") 




( ) 


( ) 




I Kill b ) 


t^Klll b j 




G45.44+0.07 


19 14 


18.291 


11 


08 


58.97 


45 


.44501 


0.06868 


0.68 


50.0 


58.5 




G45.47+0.13 


19 14 


07.362 


11 


12 


15.98 


45 


.47272 


0.13366 


8.4 


65.7 


61.9 


M 


G45.47+0.05 


19 14 


24.147 


11 


09 


43.43 


45 


.46706 


0.05321 


5.7 


56.0 


59.7 


M 


G45.49+0.13 


19 14 


11.357 


11 


13 


06.41 


45 


.49271 


0.12570 


9.7 


57.3 


60.5 


L/A 


G45.81-0.36 


19 16 


31.081 


11 


16 


12.01 


45 


.80404 


-0.35574 


14.3 


59.9 


60.3 


C 


G46.12+0.38 


19 14 


25.520 


11 


53 


Zb.yy 


46 


.11477 


0.38670 


0.90 


59.0 


55.0 


C 


G48.90-0.27 


19 22 


10.330 


14 


02 


43.51 


48 


.90213 


-0.27331 


0.52 


72.0 


68.4 


L/A 


G48.99-0.30 


19 22 


26.134 


14 


06 


39.78 


48 


.99012 


-0.29871 


0.74 


71.5 


67.5 




G49.27+0.31 


19 20 


44.859 


14 


38 


26.91 


49 


.26513 


0.31130 


8.2 


-4.5 


3.4 


M 


G49.35+0.41 


19 20 


32.449 


14 


45 


45.44 


49 


.34917 


0.41268 


11.0 


68.0 


65.5 


S 


G49.41+0.33 


19 20 


59.211 


14 


46 


49.66 


49 


.41558 


0.32588 


10.1 


-12.1 


-21.3 


C 


G49.42+0.32 


19 20 


59.817 


14 


46 


49.09 


49 


.41659 


0.32365 


3.1 


-26.0 


-21.3 


L/A 


G49.47-0.37 


19 23 


37.897 


14 


29 


59.26 


49 


.46961 


-0.37061 


2.9 


64.1 


62.4 




G49.48-0.40 


19 23 


46.191 


14 


29 


47.06 


49 


.48244 


-0.40170 


2.8 


50.0 


58.2 




G49.49-0.37 


19 23 


39.824 


14 


31 


04.94 


49 


.48936 


-0.36884 


18.9 


56.4 


60.6 




G49.49-0.39 


19 23 


43.949 


14 


30 


34.45 


49 


.48977 


-0.38750 


247 


59.2 


56.6 




G49.60-0.25 


19 23 


26.611 


14 


40 


16.99 


49 


.59933 


-0.24941 


48.6 


63.1 


57.2 


C 


G50.78+0.15 


19 24 


17.411 


15 


54 


01.60 


50 


.77881 


0.15196 


4.7 


49.1 


42.2 


s 


G52.92+0.41 


19 27 


34.960 


17 


54 


38.14 


52 


.92177 


0.41404 


4.9 


39.1 


44.8 


L/A 


G53.14+0.07 


19 29 


17.581 


17 


56 


23.21 


53 


.14183 


0.07062 


0.73 


24.6 


21.8 


L/A 


G53.62+0.04 


19 30 


23.016 


18 


20 


26.68 


53 


.61793 


0.03553 


11.2 


19.0 


23.3 


L/A 



Note. — The columns show the source name as defined in Paper I (with a prefix 'G' added to indicate derivation from 
Galactic coordinates), equatorial and Galactic coordinates, peak flux density, velocity of peak emission, the systemic 
velocity and the overall emission morphology. The positions are derived from the channel with peak emission and serve 
as the reference coordinate for the spot maps of Fig.[TJ The morphology classes are linear/arched (L/A), compact/simple 
(S), multiple (M) and complex (C). 

'^Peak flux density is 9.4 Jy after Hanning smoothing. 

'^Peak flux density is 40.8 Jy after Hanning smoothing. 
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Table 4. Absolute positions of the AMGPS 6.7 GHz methanol masers determined using 

EVLA. 



Source name RA (J2000) Dec. (J2000) I b Sp Vp V^ys Morphology 





(h m s 


) 




") 






(°) 


n 


(Jy) 


(km s ^) 


(km s 


G35.40+0.03 


18 


55 


50, 


.779 


02 


12 


19 


.08 


35, 


.39697 


0.02538 


0.66 


89.2 


95.0 


G35.59+0.06 


18 


56 


04, 


.219 


02 


23 


28, 


.34 


35, 


.58792 


0.06043 


0.20 


44.1 


49.1 


G36.02-0.20 


18 


57 


45, 


.868 


02 


39 


05, 


.67 


36, 


.01273 


-0.19719 


0.06 


92.6 


87.3 


G36.84-0.02 


18 


58 


39, 


.214 


03 


28 


00, 


.89 


36, 


.83944 


-0.02224 


2.1 


61.7 


59.2 


G37.38^0.09 


18 


59 


51, 


.586 


03 


55 


18 


.02 


37, 


.38145 


-0.08208 


0.06 


70.7 


57.7 


G37.55+0.19 


18 


59 


09, 


.985 


04 


12 


15 


.54 


37, 


.55382 


0.20089 


4.5 


83.7 


85.3 


G37.74-0.12 


19 


00 


36, 


.841 


04 


13 


19 


.98 


37, 


.73474 


-0.11193 


0.89 


50.3 


45.5 


G37.76-0.19 


19 


00 


55, 


.421 


04 


12 


12 


.56 


37, 


.75340 


-0.18915 


2.1 


54.6 


59.2 


G37.77-0.22 


19 


01 
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Note. — The columns show the source name as defined in Paper I (with a prefix 'G' added to indicate derivation 
from Galactic coordinates), equatorial and Galactic coordinates, peak flux density, velocity of peak emission and the 
systemic velocity. The positions are derived from the channel with peak emission and serve as the reference coordinate 
for the spot maps of Fig. [51 The morphology classes are linear/arched (L/A), compact/simple (S), multiple (M) and 
complex (C). 
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Table 5. MIPS 24 /xm counterparts of the AMGPS 6.7 GHz methanol masers. 
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Note. — The columns show the source name 6.7 GHz methanol maser and the equatorial coor- 
dinates of the MIPS 24 /xm point source counterpart. Only sources for which a reliable counterpart 
could be identified are shown here. Sources that are completely saturated and blanked, or which lie 
in regions of extended emission are omitted. 
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Table 6. Suggested coordinates for 6.7 GHz methanol masers that were not detected using 

MERLIN and EVLA. 



Source name 


RA (J2000) 


Dec. (J2000) 
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Note. — Based on the properties of 24 /im 
counterparts of 6.7 GHz methanol masers, the 
positions quoted here are expected to be accu- 
rate to ^ 1". 

second 24 point source is located within 
18" of the single dish position at 19''00™07^9, 
03°20'36"; see Sect. 
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Fig. 1. — Morphologies of maser spots as seen by MERLIN. For each source, the upper 
panel shows the net spectrum, while the lower panel shows the locations of the maser spots. 
Typical uncertainties in the spot positions are shown in the lower left when comparable to 
the spot symbols. The velocities of the spots are indicated adjacent to the spots themselves. 
When there is a hnear velocity trend across a set of spots, only the initial and final velocities 
are indicated at the ends of the feature. Only sources which have at least 4 maser spots 
detected at the 6a level are included in this figure. 



-33- 



G37.02-0.03 

6 
4 
2 


70 75 80 85 

LSR velocity (km/s) 



90 



6 

_ 4 
to 2 



" G37.04-0.04 

/ 


! 




J 













70 



75 80 85 90 

LSR velocity (km/s) 



95 



200 100 



AX 



-100 -200 



AX (AU) 
600 400 200 



-200 -400 



0.04 



0.02 



0.00 



-0.02 



-0.04 



77.7 o -,. . 

78.1 o 79.4 

9' 

78.2 - 79.3 



H200 



^100 



^0 



^-100 



-200 



0.04 0.02 


0.00 
Aa (") 


-0.02 -0.04 


' G37.47-0.11 1 








h 1 ' 




- |y 






55 


60 


65 70 



0.05 



-0.00 



< -0.05 



-0.10 



-0.15 



-200 



80.5,81.0 84.4 

o 847 - 85.8[»© -0 

B0.9-o, 080.3 - „ 

80.20 =80 5 %!f3:-200< 
f§l79.1 -80.3 ^ 



o 



o 



81.2 
81.8 



-400 



-600 



0.10 0.05 0.00 
Aa O 



-0.05 




LSR velocity (km/s) 



50 55 
LSR velocity (km/s) 



60 65 



1000 



0.05 



0.00 



-0.05 



AX (AU) 
500 



1.5-61.8 
63.1 



-500 



%6a8.617 

59.9 0--; 61. 5 

60.10 ^ 59.4 
59.2°o 59.5 

58.2 - 59.0 



8 

58.9 - 59.7 on 

y> <jj I 54.5 - 54.8 

a ^ 56.7-57.5 

58.0 - 58.5^ 



500 



3 

i 



-500 



1.0 



0.5 



0.0 



-0.5 



-S.O'IO^ -1.0-10* -1.5-10* 



O 49.4 - 50.3 



0.10 0.05 0.00 

Aa (") 



-0.05 



0.0 



T 10000 



-5000 



52.4 - 52.7 

S 

50.3o 



13 
< 



-5000 



-0.5 -1.0 

Aa O 



-1.5 



Fig. 1.— Contd. 
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Fig. 1.— Contd. 
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Fig. 2. — Morphologies of maser spots as seen by EVLA. As in Fig. [H the upper panel shows 
the net spectrum of each source, while the lower panel shows the locations of the maser spots. 
Typical uncertainties in the spot positions are shown in the lower left when comparable to 
the spot symbols. The velocities of the spots are indicated adjacent to the spots themselves. 
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Fig. 3. — Spectra of sources detected in the MERLIN data that have fewer than 4 maser 
spots detected at the 6a level. 
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Fig. 4. — Spectra of sources detected in the EVLA data for which maser spot morphologies 
cannot be rehably determined due to constraints of angular resolution and/or signal to noise 
ratio. 
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Fig. 5. — Histogram of the spatial extent of maser emission along the major axis of the spot 
distribution. 
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Fig. 6. — The spatial major axis extent of maser emission as a function of the velocity width 
of the maser, defined from the minimum and maximum velocities of maser emission. The 
dashed line shows a linear fit to the data after excluding outliers. 
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Fig. 7. — The spatial major axis extent of maser emission as a function of maser luminosity. 
The dashed line shows a power-law fit to the data with a power law exponent of 0.51. 
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Fig. 8. — Histogram of the separation between 6.7 GHz methanol masers and MIPS 24 /im 
point source counterparts. 



